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We study the chiral separation and chiral magnetic effects in a slab of Dirac semimetal of finite 
thickness, placed in a constant magnetic field perpendicular to its surfaces. We utilize the Bo- 
golyubov boundary conditions with a large Dirac mass (band gap) outside the slab. We find that, 
in a finite thickness slab, the axial current density is induced by helicity-correlated standing waves 
and, as a consequence, is quantized. The quantization is seen in its stepped-shape dependence on 
the fermion chemical potential and a sawtooth-shape dependence on the thickness of the slab. In 
contrast to a naive expectation, there is no chiral charge accumulation anywhere in the bulk or 
at the boundaries of the semimetal. In the same slab geometry, we also find that a nonzero chiral 
chemical potential induces no electric current, as might have been expected from the chiral magnetic 
effect. We argue that this outcome is natural and points to the truly non-static nature of the latter. 
By taking into account a nonzero electric field of a double layer near the boundaries of the slab, we 
find that the low-energy modes under consideration satisfy the continuity equation for axial current 
density without the anomalous term. 

PACS numbers: 71.70.Di, 11.40.-q, 03.65.Pm 


I. INTRODUCTION 

Nowadays there is significant interest in relativistic matter in a strong magnetic field. Such matter is intensively 
studied both experimentally and theoretically. It has a number of applications in high-energy physics and astrophysics 
(e.g., in the context of compact stars, heavy-ion collisions, and the early universe), as well as in condensed matter 
physics (e.g., in the context of novel Dirac/Weyl materials). 

In condensed matter physics, after the first reports of three-dimensional Dirac semimetals appeared two wars ago 
[in, the field exploded with numerous investigations of their exotic properties (for reviews, see Refs. (JQ)- 
of the key aspects of such materials is a well-defined chirality of the low-energy quasiparticles, described by Weyl 
fermions. Because of the chiral anomaly Q, the chirality is not a conserved charge and, thus, the corresponding 
quasiparticles must come in both chiralities In Dirac semimetals, the low-energy spectra of the quasiparticles 
of opposite chiralities are degenerate. Such degeneracy is often protected by symmetries (e.g., the symmetry under 
time-reversal or parity). If the corresponding symmetry is broken, however, a Dirac semimetal may turn into a Weyl 
semimetal, in which the degeneracy of the quasiparticles with opposite chiralities is lifted. A number of materials of 
this latter type has been recently reported as well [mi. 

Phenomenologically, a chiral asymmetry in relativistic matter may be introduced via a nonzero chiral chemical 
potential [l^. Such a chemical potential couples to a difference between the number densities of the left- and right- 
handed fermions and enters the Lagrangian density through the term where 'if = The latter produces 

a chiral asymmetry in magnetized relativistic matter and leads to a nondissipative electric current j = e^B/i 5 /( 27 r^) 
in the presence of an external magnetic field B flfij - fisjl (see also Refs. fTol-l^). This phenomenon is known in the 
literature as the chiral magnetic effect (CME) and its origin is related to the famous chiral anomaly Q- Moreover, 
the charge-dependent correlations and flow, observed in heavy-ion collisions at RHIC [23 - [^ and LHC [^, appe ar 
to be in a qualitative agreement with the predictions of the CME [2^, (for recent reviews, see also Refs. |3ll - [^ ). 
In the context of condensed matter physics, it was also suggested that the measured quadratic field dependence of 
the magnetoconductance in ZrTes is an indication of the chiral magnetic effect [s^ . 

Unlike the chiral chemical potential, which is a rather exotic quantity and is not very well defined theoretically, 
the chemical potential /i (associated, for example, with a conserved electric or baryon charge) is common in many 
physical systems. It was shown in Refs. |35M37j | that a nondissipative axial current density js = —eB/r/(27r^) exists 
in the equilibrium state of noninteracting massless fermion matter in a magnetic field. This effect is known as the 
chiral separation effect (CSE). In fact, as suggested in Refs. [3^ the CSE may lead to a chiral charge separation 
(i.e., effectively inducing a nonzero chiral chemical potential fj-b) and, thus, trigger the CME even in the absence of 
topological fluctuations in the initial state. 
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The physical and mathematical reasons for the chiral asymmetry in relativistic matter in a magnetic field are quite 
transparent (for an elegant exposition, see also Ref. [l^). In essence, its origin is connected with the spin-polarized 
nature of the lowest Landau level (LLL). The corresponding fermionic states are also characterized by a well-defined 
longitudinal momentum and, thus, chirality. Moreover, the states with opposite signs of the longitudinal momenta 
carry opposite chiralities and, thus, lead to a nondissipative axial current density js = —eB^/(27r^) [H, H^. 

It was argued in Refs. that nondissipative currents in magnetized relativistic matter are determined by 

the topological currents induced exclusively in the LLL and are intimately connected with the chiral anomaly. This 
fact is directly connected with the well-known result that the chiral anomaly in a magnetic field is also generated 
exclusively by the LLL [4l|. The first studies of interaction effects on the chiral asymmetry of relativistic matter in 
a magnetic field were performed in Refs, [s^ by using Nambu-Jona-Lasinio models with local interaction. In 

particular, it was found that the interaction unavoidably generates a chiral shift A [sM when the fermion density 
is nonzero. It enters the effective Lagrangian density through the following quadratic term: when the 

magnetic field is directed along the z direction. The meaning of the chiral shift parameter is most transparent in the 
chiral limit: it determines a relative shift of the longitudinal momenta in the dispersion relations of opposite chirality 
fermions, —>■ ± A, where the momentum is directed along the magnetic field. Furthermore, as shown in 
Refs. [ 33 , E3jEJ|i th® chiral shift A is responsible for an additional contribution to the axial current density. Also, 
such a dynamically generated chiral shift splits each Dirac node into a pair of Weyl nodes of opposite chirality, thus 
producing a Weyl semimetal from a Dirac one [d^. Recently, another interesting mechanism for inducing a nonzero 
chiral shift was proposed Eli E3l • It uses a circularly polarized light and works even at zero chemical potential. 

Usually the chiral magnetic and chiral separation effects are considered in the literature in unbounded material 
media. In practice, however, all physical systems (except for the CME and CSE in the early universe, perhaps) are 
finite. It is natural then to ask about the role of boundaries and finite-size effects in the chiral magnetic and separation 
effects. Indeed, even if one assumes that, in the bulk of a bounded medium, the electric and axial currents are the same 
as those in infinite systems, they should get modified near the boundaries. This should be an immediate consequence 
of the continuity equations if the currents are required to vanish outside the material. This simple observation was 
the main motivation for the present work. Here we will study the chiral separation and chiral magnetic effects in 
a slab with an external magnetic field perpendicular to the boundary planes of the slab, which yields the simplest 
realization of a finite-size system. 

The paper is organized as follows. In Sec. [Ill we introduce the model of a slab with the Bogolyubov boundary 
conditions [i^ . when the mass (band gap) parameter in the vacuum outside the slab is taken to be the largest mass 
(energy) parameter in the model. Note that we will use both the terms “mass” and “band gap” interchangeably 
throughout the paper. Section Hill is devoted to the analysis of the axial current density (i.e., the CSE) in the slab 
with such boundary conditions. The CME is considered in Sec.|TVl In Sec.|V]we discuss the chiral anomaly in vacuum 
regions near the surface of the slab. The discussion of the main results is given in Sec. IVII Some technical details, 
including the derivation of the Landau-level wave functions and the implementation of the boundary conditions, are 
presented in Appendices |A] and [B] respectively. 

Throughout this paper, we set h = 1 and c = 1. 


II. MODEL 

The Hamiltonian of the low-energy model of a Dirac semimetal slab situated between the planes z = —a and z = a 
reads 


H = 


d^r d>^(r) [vpa. ■ (—iV + eA) + 7 °m(z)] 'l'(r), 


where a = 7 ^ 7 , and 7 are the Dirac matrices in the chiral representation, i.e.. 


7 


0 




7 = 


0 tr] 

-cr 0 J ’ 


7 ® = i7°7^7^7^ 



( 1 ) 

( 2 ) 


Here I2 is the two-dimensional unit matrix and cr = (cri, (72, <73) are the Pauli matrices. By assuming that the external 
magnetic field B is directed along the z axis, we will find it convenient to use the vector potential in the Landau 
gauge, i.e., A = (0, i?a:,0). The other notations are as follows: e is the electron charge, vp is the Eermi velocity, 
and m(z) = M 9[z‘^ — a?) + m9{a? — z^) is the Dirac mass (band gap) function [here 9{x) is the unit step function]. 
The case of a Dirac semimetal with a zero band gap in the bulk is obtained by taking the limit m —>■ 0. In the 
model at hand, we assume that the “vacuum” gap parameter M is much larger than all characteristic energy scales 
in the slab (e.g., the work function and/or relevant quasiparticle energies). Interestingly, such a model of the slab is 
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nothing else but a generalized Bogolyubov bag model [d^. In the studies of graphene, similar boundary conditions 
with an infinite gap outside the material are known as the infinite mass boundary conditions [50l | . The same 

idea, albeit with a finite-size band gap, is also utilized in modeling a potential barrier in the context of the Klein 
paradox in graphene; see Ref. [5l|. While m = 0 was used in the original Bogolyubov model, we will treat to as a 
free parameter in the analysis below. (Note that the m 0 case is of interest not only from a theoretical viewpoint, 
but could be also investigated experimentally in a more general class of Dirac semimetals/metals, e.g., such as a 
bismuth alloy Bii_ 2 ;Sba; at small concentrations of antimony [^. [s^. where the Dirac gap is nonzero.) For a good 
review of bag models in hadron physics, see Ref. [l^. Generically, in all such models, the hadrons are described as 
bags with massless fermions (quarks) confined inside. In order to prevent massless fermions from leaving the bag, 
one requires that the normal components of the momenta and hence the currents across the surface vanish. Such 
boundary conditions necessarily break the chiral symmetry for fermions [s^. From the physics viewpoint, this is 
unavoidable because massless quasiparticles experience a helicity flip (and, thus, a chirality change) whenever the 
directions of their momenta change due to scattering from the boundary. This fact will be crucial for our analysis 
below. In particular, in the Bogolyubov model with the vanishing gap in the bulk, to —>■ 0, the chiral symmetry is 
explicitly broken by the inclusion of the (infinitely) large vacuum gap parameter M. 


III. CHIRAL SEPARATION EFFECT 


In this section, we calculate the axial current density and study the CSE in a Dirac semimetal with a slab geometry. 
By making use of Hamiltonian o, the ground state of the system will be obtained by filling all quasiparticle states 
with the energies less than the Fermi energy Ep = where /r is the chemical potential. 

Let us start by determining the energy spectrum and the electron wave functions for a Dirac semimetal slab in 
a constant magnetic field. Before proceeding to the slab case, however, it is instructive to start from presenting 
the Landau-level wave functions in an infinite space. By making use of the chiral representation Q, we derive the 
following results for the wave functions (see Appendix for details): 
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where I = l/y/\eB\ is the magnetic length, ep = is the Landau energy scale, and ^ = x/l+pyl. For convenience, 

here, we fixed the sign of electric charge so that sign(ei?) = -1-1. Additionally, we introduced the following harmonic 


oscillator wave functions: K„(0 = 


p-C^/2 


where iLn(0 ^'^e the Hermite polynomials. (For the Landau-level 


wave functions in the standard representation of the Dirac matrices, see Ref. [55j.l The corresponding Landau-level 
energies are given by En = ± VpP^ + m? -f 2ne|. 

Now, in the case of a slab with a finite extent in the z direction, for every plane wave with a wave vectoror pz, 
propagating in the positive z direction, there should be also a plane wave with a wave vector —pz, propagating in the 
opposite direction. Therefore, the general solution for the nth-Landau-level wave function in the slab is given by a 



























4 


superposition of two counterpropagating plane waves, or equivalently, a standing wave: 
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These wave functions inside the slab should be matched to the corresponding solutions in the vacuum. This is done 
in Appendix |B] using the Bogolyubov bag model, in which wave functions outside the semimetal satisfy the Dirac 
equation with an infinitely large vacuum gap parameter M. By enforcing the boundary conditions, we also find that 
the wave vector should satisfy the following spectral equation: 


vfPz cos {2apz) + rn sin {2apz) =0, (7) 

where Pz ^ 0. The final expressions for the Landau-level wave functions in the slab are given in Eqs. dni), dug), 
and (IB16I) . In essence, they have the form of standing waves with discrete wave vectors Pz that satisfy Eq. ([7]). It is 
important to note that, while there are two independent solutions for each wave vector in the higher Landau levels, 
there is only one independent solution for each wave vector in the LLL. 


A. Axial current density 


In this subsection, we calculate the axial current density inside the slab. In terms of the Landau-level wave functions, 
the corresponding ground state expectation value is given by 


01) = f f/(PO^^f’^'IlabO)n=07'’7V^slab(r)„=0 + Xl£/OO^^f’^slib(l')n7°7V«'l)ibO)") > 

Pz \ i=l n=l / 


( 8 ) 


where the contributions of quasiparticles from both valence and conduction bands are taken into account via the use 
of the following generalized distribution function: 


f{Pz) = 


1 


I 


(y7fpf+m2+2n^-/i)/T ^ ^ /T 


( 9 ) 


I 


Here p is the chemical potential (Fermi energy) measured from the Dirac point. This distribution function accounts for 
the fact that the quasiparticles of the conduction and valence bands carry opposite charges. In the zero-temperature 
limit, which we use in the following, this function simplifies: 


hm f{pz) = sign(^)6» - vj,pl - - 2nel) , 

T—^0 


( 10 ) 


where 0(x) is the unit step function. 

When calculating the axial current density, it is instructive to separate the contribution of the spin-polarized lowest 
Landau level from the contributions of the higher Landau levels (n > 0). The LLL contribution is obtained by making 
use of the wave function in Eq. (Eg. The zero-temperature result reads 
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where we also made use of the spectral equation in Eq. ©■ As we see, in a general case when to 7 ^ 0, the LLL 
contribution to the axial current density has a nontrivial dependence on the z coordinate. Here it may be appropriate 
to mention that the axial current density is well defined even in the gapped case when the axial charge is not conserved. 
In fact, even in the extreme non-relativistic limit, it has a clear physical meaning as a spin polarization. 

As expected, in the chiral (gapless) limit, the axial current density is independent of the z coordinate and the 
explicit result reads 




0,m^0 — 


eBvp sign(/i) 
doTT 


2 2 2 \ eBvFsign{^i) 

2^0 [n - vpp^) = -^- K 

Pz 


4a7r 


( 12 ) 


where we took into account the spectral equation which reduces down to cos {2apz) = 0 in the gapless case. The 
latter also implies that the allowed values of the wave vector are p^^\ = {2k — l) 7 r/( 4 a), where fc is a positive integer. 
Because of the unit step function in Eq. the result of the sum is given by fcmax = [‘2a\P'\/{vpTr) + 1 / 2 ] where [...] 
represents the integer part. As expected, in the limit of a —?► 00 , the above result reduces to the well-known relation 
for the chiral separation effect in an infinite system, i.e., (j|) = —\eB\p/{2Tr^). 

The result in Eq. m shows that, in a slab of finite thickness a, the axial current density is quantized. This is a 
qualitatively new feature that did not exist in an infinite system. It is a natural outcome of having the wave functions 
in the form of standing waves and the quantization of the wave vector in a slab geometry. As is easy to see from 
Eq. (I12L the height of the steps in the axial current is proportional to the magnetic field and inversely proportional 
to the thickness of the slab, i.e., S{j^) = eBvp/{4:aTT). When (j|) is plotted as a function of ap/vp, the widths of the 
steps are given by tt/2. It would be interesting to explore whether such a quantization can be observed in experiment. 

In the case of higher Landau levels (n > 0), there are separate contributions to the axial current density from each 
of the two independent modes in the slab, i.e., 'I'g]',]|^(r)„ and 'I'gj^,]|^(r)„. By making use of the explicit expressions 
for the corresponding wave functions in Eqs. (IB15I) and (IB16|) . we find that, for each p^, the two contributions have 
opposite signs, i.e.. 




eB vpp^ {m? + Vppl) [1 - cos (2p^z) cos {2p^a)] 
27r 2apzE^ + vppzm — sin {4apz) 
eB vppz (to^ -I- Vppl) [1 - cos {2pzZ) cos (2pza)] 
2 tt 2apzE^ + vppziTi — nej^ sin (dap^) 


(13) 

(Id) 


Therefore, there is no net contribution to the axial current density due to higher Landau levels. In other words, just 
like in the case of an infinite space , the axial current density in a semimetal slab is determined exclusively by the 
LLL contribution (|TT|l . 

Our numerical results for the axial current density (fTT|l in a slab geometry are presented in Eigs. [U O and m For 
comparison, in Eigs. [1] and [2] we also show the results for an infinite space. The dependence of the dimensionless axial 
current density 2tt‘^{ jf)aP/vp, in the middle of the slab (z = 0 ), as a function of the chemical potential is presented 
in Eig. [H [Recall that the axial current density in Eq. (ED is not uniform in space when to 7^ 0.] The three panels 
show the results for three representative values of the band gap: am/vp = 0 (left panel), am/vp = 2 (middle panel), 
and am/vp = 6 (right panel). To plot the figure, we used uf/o = 25 meV. As expected, a nonzero current density is 
obtained only when p > to, i.e., when the chemical potential is larger than the size of the band gap. In agreement 
with the earlier observation, we also see that the axial current density is quantized in the slab. 



FIG. 1; (Color online) The dimensionless axial current density in an infinite space (red dashed line) and in the middle of the 
slab (blue solid line) plotted as a function of chemical potential for the three values of the band gap; amjvF = 0 (left panel), 
am/vp = 2 (middle panel), and am/vp = 6 (right panel). To plot the figure, we fixed vp/a = 25 meV. 
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FIG. 2: (Color online) The dimensionless axial current density in an infinite space (red dashed line) and in the middle of the 
slab (blue solid line) plotted as a function of the width of the slab a for the three values of the band gap: m//r = 0 (left panel), 
mill — 0.25 (middle panel), and mlii = 0.75 (right panel). To plot the figure, we fixed vpjii— 12.5 A. 


The dependence of the dimensionless axial current density 2'k‘^ jii, in the middle of the slab (z = 0), as a 
function of the width is shown in Fig. [21 The three panels show the results for the following three values of the band 
gap: m/ll = 0 (left panel), m/ii = 0.25 (middle panel), and m/ii = 0.75 (right panel). To plot the figure, we used 
Vf/ fJ- = 12.5 A. As we find, the functional dependence of the current density on the width of the slab has a sawtooth 
shape. This is rather natural consequence of the quantization of the wave vector. When the size a becomes large, 
the finite-size effects quickly decrease and the result approaches the limit of an infinite system. We also see that a 
nonzero gap has a damping effect in the dependence on a. [Notice the difference in the vertical scales in the three 
panels of Fig. |2j] 

It is interesting to note that, in the gapped case (m ^ 0), the steps in the axial current density have different 
heights, see Fig. |TJ This is pronounced the most in the lowest few steps of the current density. (Note, at the same 
time, that the widths of the steps in aii/vp remain nearly, although not exactly, the same.) The corresponding steps 
are determined by the low-energy quasiparticle states with the smallest few wave vectors, i.e., the wave vectors which 
are modified the most by a nonzero gap. This can be checked explicitly in the limit of a small but nonzero band gap. 
In such a limit, an analytical expression for the wave vectors can be obtained approximately by solving the spectral 
equation 0. The result reads ~ (2fc — l)7r/(4a) -I- 2 m/[7r?;_F(2fc — 1)] where fe is a positive integer. By making 
use of this, one can check that the contributions from the states with small k vary a lot because of great variations in 
the values of cos(2p^™^a); see Eq. (ITT]) . (Away from the middle point in the slab, the heights of the steps also vary. 
This is clear from Fig. jT]) With increasing the value of a\ii\/vF, on the other hand, the effect of the gap diminishes 
and the functional dependence of the axial current density gradually approaches the result in the gapless limit. This 
is understandable since the states with large wave vectors, which dominate the result in such a regime, are insensitive 
to the size of the gap. 



FIG. 3: (Golor online) The dimensionless axial current density plotted as a function of the chemical potential and z/a for the 
three values of the band gap: am/vp = 0 (left panel), am/vp = 2 (middle panel), and am/vp = 6 (right panel). To plot the 
figure, we fixed vp /a = 25 meV. 

As is clear from Eq. in the case of a nonzero band gap, the axial current density is not uniform: it depends 
on the position 2 inside the semimetal. The corresponding dependence is shown explicitly in Fig. [31 where we present 
the results for the gapless case am/vp = 0 (left panel) alongside with the results for the two cases with nonzero band 
gaps: am/vp = 2 (middle panel) and am/vp = 6 (right panel). The fact that the axial current density is constant 
inside the slab in the gapless case is consistent with the continuity equation. Indeed, in the case of the vanishing 


















7 


gap and in the absence of electric fields, the axial current is conserved inside the slab. In contrast, there are clearly 
visible oscillations of the axial current density in the other two cases when m ^ 0 (see the middle and right panels 
in Fig. [3]). Taking into account that the axial current is not conserved in the gapped case, such oscillations are not 
forbidden. Moreover, as we see, the larger is the gap, the more pronounced are the oscillations. When the gap is 
fixed, we also find that the local amplitude of the axial current density oscillations is not the same across the whole 
sample: it decreases with the increasing distance from the slab boundaries. 

It may be interesting to emphasize that, in a finite-thickness slab, the axial current density originates from helicity- 
correlated standing waves. Of course, the underlying roots of this are (i) finite thickness of the slab and (ii) the 
spin-polarized nature of the LLL. This can be seen explicitly from the structure of the LLL wave function in Eq. (iBTl) . 
which takes the following simple form in the gapless limit: 


^slab(l’)n—0 — 


^iPyV-{lPy+x/l)'^/2 



(15) 


Such a standing wave is made of a pair of counterpropagating plane waves carrying opposite chiralities or helicities. 
This is clear since, in the chiral representation used, see Eq. (ED, the upper (lower) two components of the wave 
function describe right-handed (left-handed) modes. The configuration in Eq. dTS]) can be interpreted as follows: the 
helicity of each plane wave propagating to the boundary flips the sign after the reflection. This is indeed expected 
for the spin-polarized LLL states. A nontrivial helicity correlation of the LLL standing waves remains even in the 
gapped case, but it is not as transparent. Such a correlation is the key feature that is responsible for the axial current 
density in a semimetal slab. 

Before concluding this subsection, let us briefly discuss the implications of the chiral separation effect in a slab of 
finite thickness. Naively, one may expect that the axial current density in the bulk of the semimetal should lead to 
an accumulation of positive chiral charge on one side of the slab, and negative chiral charge on the other. As is easy 
to check, however, this does not occur. In fact, the axial charge density is identically zero everywhere, = 0. 


B. Axial current in the vacuum with finite band gap 


As the results in the previous subsection show, in a slab geometry, there is a nonzero axial current density in 
the bulk of semimetal. However, the current should be vanishing in the vacuum because of the imposed boundary 
conditions. Then, the current density should presumably go to zero in the surface layer. In order to clarify this issue, 
in this subsection, we investigate the continuity equation for the axial current density at the vacuum side of the slab. 

In the absence of an electric field, the axial charge and current densities should satisfy the following continuity 
equation Q: 




(16) 


in all regions of space, including at the boundary of the slab. Note that, in order to have the continuity equation 
well-defined, we will assume that the vacuum band gap M is large, but finite. The corresponding results for the LLL 
wave functions are derived in Appendix IIII B[ see Eqs. (IB22I) . (IB23I) . and (IB24I) . Also, in the same appendix, we 
derive a modified version of the spectral equation for the wave vector see Eq. (IB21I) . 

By making use of the wave functions in Eqs. (IB23I) and (IB24|) . we can also compute the axial current density outside 
the slab. In the region z > a, the corresponding result is given by 
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where the sum is performed over the discrete values of the wave vectors that satisfy Eq. (IB21|) . The result in the 
region z < —a is similar, but one should replace —{z — a) —>■ (z -I- a). 

In order to check the continuity equation in the regions outside the slab, let us calculate the ground-state expectation 



























value of the quantity that appears on right-hand side of Eq. (1161) . It is straightforward to show that 
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In the region z < —a, the corresponding quantity is obtained by replacing —{z — a) —> (z -|- a) and flipping the overall 
sign. 

By taking the derivative of Eq. (HD with respect to z and making use of Eq. dUD, we check that the continuity 
equation (ITSl) is indeed satisfied in the regions outside the slab. Therefore, the axial current flows to the boundary, 
where, due to the presence of the large vacuum mass M, it exponentially vanishes. Moreover, as claimed earlier, this 
happens without any chiral charge accumulation at the boundaries. 


IV. CHIRAL MAGNETIC EFFECT IN THE SLAB 


In this section, we will study the chiral magnetic effect in a slab geometry (for simplicity, only the case of the chiral 
limit will be considered here). In order to do this, we need to introduce a nonzero chiral chemical potential /rs into 
our model of a semimetal. In the chiral limit, one might try to introduce /rs by just replacing the fermion number 
chemical potential /r with in the distribution function (1101) for the right- and left-handed particles, respectively. 

By recalling, however, that in the presence of the boundaries the chiral symmetry is broken even in the case of gapless 
Dirac fermions in the bulk, we know that does not correspond to a conserved quantity. In such a situation, it is 
more convenient to treat as a phenomenological parameter that modifies the model Hamiltonian, i.e., 

H = f (Pr \I'(^(r) [vpa ■ (— iV -f eA) -/r - M57"]'l'(r). (19) 


It is worth noting, that in the model under consideration /rs determines the relative energy shift of the left- and 
right-handed Weyl nodes. We calculate the ground-state expectation values by using the Schwinger prescription, 
where the summation over energy eigenvalues is performed with the distribution function fo{E) = —sign(E)/2. [Note 
that, in this case, the effective distribution function /o(i?_) -|- /□(£'+), where E± = —fi ± y/vppl + gives the 

same distribution as Eq. (ITUl) in the gapless limit.] 

For Hamiltonian (|19L we find the following LLL wave function inside the semimetal: 



' 

/0\ 




^siab(r)„=o = YoiOP^yy < 

< 

O O 1-^ 

+ PP^^Co 

1 

\l) 

] 


( 20 ) 


where Pz = —{E + p pP)Ivf and pz = {E + p — pp)/vF- The corresponding wave functions outside the slab are 
given by Eqs. m and (IB2|) . By making use of the matching conditions in Eqs. (IB3I) and (IB4I) . we derive the spectral 
equation for the wave vectors, cos[a(p^ ~ Pz)] = 0. By solving this equation, we find that the energy parameter E 
can take only the following discrete values: 


Ek 



( 21 ) 


where fc is a positive integer. After enforcing the boundary conditions, we also derive the following final expression 
for the wave function inside the slab: 


0 

giPzZ 

0 

and find that the corresponding LLL contribution to the electric current vanishes, 

OO 

if) ='^f{Ek)vF'^lubi^)n=0j°l^'^sUh{r)n=0 = 0 . 

k^l 


lkslab(l')n—0 — 


Yoi^PPyy 

2'/aI 



(23) 
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Although this result is natural in the Bogolyubov model, it may appear surprising because it is precisely the LLL 
contribution that saturates the CME in an infinite system. To complete the analysis of the CME in a finite slab, 
however, we should still analyze the contributions of higher Landau levels (n > 0). For a semimetal described by 
Hamiltonian (I19L the wave functions in higher Landau levels are given by 


llfslab(l*)n — C 


WyV 


^ -ivF (+ + p.) K-i(e) ^ 




0 

0 


/ 


+ 026 


iP. 


\ 


^ -IVF + P.) Yn-m ^ 


036 


-iP^z 


V^2ne|r„(0 

0 

0 


Oi6 


— iPzZ 


zvF ( \IP^ + - p,) 

\ 


^2n6\Yn{i) 

0 

0 


^VF[JP^+^-^+Pz]Yr..3{0 


V 


v'2ne|r„(,f) 


/ 


(24) 


where P^ = yJ{E + ^ + — 2ne\ and = Vp^ yJ{E + ^ — fi^)^ — 2ne\. 

The corresponding spectral equation is obtained by matching the wave functions in the bulk with those in the 
vacuum; see Eqs. (IB3I) and (lB4l) . Its explicit form reads 


(l_g4*aP.) [(p + ^)2 _ ^2 _ 2 „g 2 j j'^^g4*aP. ^g4zaP. ^ 4 g 2 *a(P.+A) ^g4*a(P.+A)^ ^ q_ (25) 

When /X 5 yf 0, we could solve this equation only numerically. By making use of the spectral equation, we can write 
down the formal solution for the semimetal wave function in the slab. The corresponding expression is bulky and not 
very informative. Therefore, we do not present it here. Instead, by making use of the numerical solutions to Eq. p5|) , 
we calculate the contribution of the higher Landau levels to the electric and axial current densities. Both results 
vanish within the numerical precision used. 

Thus, unlike the axial current density of the CSE, the electric current density of the CME is absent in the slab. 
This result might have been expected from general considerations. Because of the Bogolyubov bag model boundary 
conditions (with an infinite band gap in vacuum), there should be no electric current across the boundary. Taking 
into account that the electric current is not anomalous, this means that the continuity equation for the electric current 
can be satisfied only if the electric current vanishes inside the slab. This is exactly what the direct calculations give. 
Moreover, in the model at hand, the electric current will remain vanishing even in the limit a —>■ oo. This fact is 
consistent with the case of infinite systems in equilibrium considered in Refs. [sgI - Is^ . 

In passing, let us briefly discuss the simplest case of a Weyl semimetal. The corresponding model Hamiltonian 
will be similar to that in Eq. (HH), but will include an additional term, —vf f d^r 4 'f(r) (a • A)7®4t(r). The value 
of A = Abz, which is often called the chiral shift parameter, determines the separation between the Weyl nodes 
in momentum space. Here, for simplicity, we assume that A points in the +z direction. This may be sufficient 
because the components of the chiral shift parallel to the slab surface are not expected to modify the currents in the 
slab. When A 0, we find that the LLL wave functions in the bulk have the same form as in Eq. (1201) . but with 
Pz = —{E + /i + P 5 )/vf + A and pz = {E + p — P 5 )/vf — A. Then, by making use of the wave functions outside the 
slab, see Eqs. m and (IB2I) . and enforcing the boundary conditions in Eqs. (|B3I) and (|B4I) . we hnd that the energy 
can take only the following discrete values: 


Ek = -p + vfA + “ 0 ’ 

where fc is a positive integer. Compared to the result in Eq. (ED, the only difference here is the change: p — >■ pl — vfA. 
This finding is in agreement with the result in an infinite system, obtained in Ref. . The calculation of the higher 
Landau level contributions to the axial current is more challenging in the general case ^ 0. Therefore, we restrict 
our discussion to the simpler case, = 0. In contrast to the case of an infinite system, we find that the contributions 
of higher Landau levels to the axial current density vanish. Therefore, our results for the axial and electric current 
densities in a Dirac semimetal slab, see Eqs. (fT^ and (|23)l . remain qualitatively the same also in the case of a Weyl 
semimetal slab, but with the replacement p ^ fi — vf A. 
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V. WORK FUNCTION, ELECTRIC FIELDS, AND CONTINUITY EQUATION FOR AXIAL CURRENT 


The work function of a solid is defined as the energy needed to remove an electron from the solid to vacuum (see, 
e.g., H). Microscopically, it can be thought of as the result of a “confining” electric field £ near the surface, resulting 
from the electron density leaking slightly out of the material. In the problem at hand, we consider a seinimetal slab 
with a nonzero magnetic field perpendicular to its surface. This means that £ • B ^ 0 near the surfaces. As we know, 
such a field configuration has a nontrivial (anomalous) contribution to the continuity equation for the axial current 
0 , i.e., 


doji + V • js = • B. 


(27) 


Therefore, it is important to investigate the role of this relation near the surfaces. 

Let us assume that the electric field £ is perpendicular to the surface and exists only in a thin layer of depth Xe- 
We can conveniently describe such a configuration by using a time-like component of the vector potential, Aq, i.e., 

— a — Xe < z < —a : Aq = —£z, (28) 

a < z < a + Xe '■ Aq = £z. (29) 

By assumption, Aq vanishes inside the slab (i.e., — a < z < a), as well as in the vacuum regions outside the thin 

surface layers (i.e., for z < —a — Xe and z > a + Xe)- 

For the sake of clarity and brevity, we consider the case of zero gap inside the semimetal, to = 0. Inside the 
semimetal, where electric field is absent, the LLL wave function is given by 




a<,z<.a 


= e 


— iEt 


YoiOe 


iPyV 


( \ 

0 

V / 


(30) 


Although the Dirac equation admits exact analytic solutions in constant electric and magnetic fields [b^, we find it 
more convenient and transparent to obtain the corresponding solutions in the region near the slab boundaries, where 
electric field is present, in the first order of the perturbation theory in electric field £. Since we treat electric field in 
perturbation theory and only the LLL modes contribute to the axial current in a magnetic field, we begin with the 
following ansatz [compare with the LLL wave functions (l30ll in the absence of electric field]: 


'^vac 




( ° \ 

1 

0 

V (t)2{z) ) 


(31) 


Then the Dirac equation in the region a < z < a -\- Xe implies the following equations for the functions (ji and 4i2- 


E(j)i{z) — vpiidz + e£z/vF)4'i{z) + M(j) 2 {z) = 0 , ( 32 ) 

E(f) 2 {z) + vpiidz — e£z/vF)4>2{z) + M(j)i{z) = 0 . ( 33 ) 

We seek solutions of this system of equations as the Taylor series in £ retaining only its two first terms 

(jiiz) = foiz) - e£fi{z), ( 34 ) 

h{z) = go{z) - e£gi{z). ( 35 ) 

By substituting Eq. (1551) into Eq. (155)) . we obtain the solutions 

fo{z) = Aie~P‘^^ + A 2 eP-\ g^^z) =Aie~P-^{E + ivFp'z) + - ivFp'z) , ( 36 ) 


fi{z) = ^3 {Aie-P> [i{vj-p'^ + 2M‘^z) - 2 E^z^Ivf + E{vf + 2vfpIzZ + 2[Mz)‘^Ivf) - 2iE^z] + 

\i{v^FP'z ~ ‘2M^z) + 2E^z^ /vf — E{vf — 2vfp'zZ 2{Mz)‘^/vf) + 2iE‘^z] | , ( 37 ) 

gi{z) = ~ [-2E^z'^{E + ivFp'z)/vF + M'^{vf + 2vfp'zZ + 2Ez'^{E + ivFp'z)/vF)] + 

- 1 -^ 26 ^*^ [2E^z'^{E — ivFp'z)lvF — M^{vf — 2vfp'zZ + 2Ez^{E — ivFp'z)/vF)] | , ( 38 ) 
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where p'^ = Vp^y/M'^ — E'^ and Ai and A 2 are constants. By matching the wave function (l30l) and the vacuum 
functions 


/ 


Ip-a- 


XE<z<—a 




0 

fo{z + a) + e£fi{z + a) 
0 

\ go{z + a) + 6 ^ 51 ( 2 ; + a) 

( 0 

foiz- a)- e£fi{z- a) 
0 

V goiz - a) - e£gi{z - a) 
( ° 

\ -i 


i^a<z<a+x, = e-*^‘yo(Oe" 


(39) 


(40) 


(41) 


(42) 


at z = —a — Xe-, z = —a, z = a — Xe, and z = a, and taking into account the normalization conditions, we find that 
the linear in £ terms of the wave functions dini) and (HD) do not contribute to the continuity equation for axial current 
because 


ds{jl)-a-x.<z<-a - = 0{e^£^), (43) 

d3{jl)a<z<a+x, - J = 0{e^£^). (44) 

This result is quite natural because the anomalous term in the chiral anomaly is connected with the lack of a chiral 
symmetry invariant regularization of the famous linearly divergent triangle diagram Q- Thus, it is high-energy 
modes whose contribution is divergent and should be regularized which are responsible for the anomalous term in the 
continuity equation for axial current. Consequently, we conclude that the low-energy modes confined in the semimetal 
satisfy the continuity equation for axial current without the anomalous term. 


VI. CONCLUSION 

In this paper, we studied the chiral separation and chiral magnetic effects in a Dirac semimetal with a slab geometry 
placed in a constant magnetic field perpendicular to its surfaces. We used the Bogolyubov boundary conditions |48l | 
at the surfaces of the slab. It is worth mentioning that such a model was originally used in high-energy physics for 
the description of hadrons within the framework of the bag models of quarks [54| . This model assumes that the 
quasiparticles have a large band gap (Dirac mass) in the vacuum regions outside the slab. 

Using this model setup, we derived analytically the spectral equation and the wave functions for the bulk modes. 
Furthermore, we calculated the axial current density and found that, just like in an infinite system, only the LLL modes 
contribute to this quantity. We show that the main consequence of a finite slab thickness is that the axial current 
density becomes quantized. The corresponding quantization could be revealed in its functional dependence on the 
chemical potential and the thickness of the slab. The underlying reason for the quantization is the discretization of the 
wave vectors and energy levels in the slab. In other words, the axial current density originates from helicity-correlated 
standing waves, associated with the LLL. It is also interesting to point that there is no chiral charge accumulation 
on the semimetal boundaries, as might have been naively expected. The size of the quantization steps in the axial 
current density as a function of the chemical potential, see Fig.[Tl is given by S{j^) = eiIu_F/(4a7r). Because the steps 
are inversely proportional to the thickness of the slab a, most likely the corresponding structure can be observed in 
experiment only if the samples are sufficiently thin. 

As we show in this study, the dependence of the axial current density on the thickness of the slab a has a very char¬ 
acteristic sawtooth shape, see Fig. [21 With increasing the value of a, the quantization effects become less pronounced 
and the axial current density gradually approaches the limit of an infinite system. Interestingly, the quantization 
persists even in the case of Dirac quasiparticles with nonzero gaps. However, a nonzero gap has a damping effect in 
the dependence on a. 







12 


We also find that, in the case of gapped Dirac fermions, m 7 ^ 0, the axial current density is not uniform across the 
semimetal slab, but oscillates as a function of the position, see Fig. [3] This is in contrast to the chiral limit when the 
axial current is constant. In fact, the larger is the gap, the more pronounced are the oscillations. Formally, such a 
behavior is connected with a nonconservation of the axial charge. It appears that the amplitude of the oscillations 
decreases with increasing the distance from the boundaries. To the best of our knowledge, the possibility of such 
oscillations and their key features have not been reported before. It would be very interesting to test this scenario in 
the lattice simulations. Also, this phenomenon could potentially be investigated experimentally in a semimetal with 
a spontaneously generated gap, e.g., via magnetic catalysis , or another mechanism. 

By introducing a nonzero chiral chemical potential fj,^, we also analyzed the chiral magnetic effect in a semimetal 
slab. In contrast to the axial current density due to the CSE, the corresponding electric current due to the CME is 
absent in a semimetal slab. In retrospect, the vanishing result is a natural outcome of the non-anomalous continuity 
equation for the electric current and the Bogolyubov boundary conditions [b^. Indeed, the latter are equivalent to 
the requirement of vanishing electric current from the semimetal to vacuum. It is worth emphasizing, however, that 
in this study we limited our consideration only to static configurations. This automatically excludes all transient 
phenomena in which nonzero electric currents are generated out of equilibrium and evolve in time. The corresponding 
generation of the CME in finite-size samples would be very interesting and should be studied in the future. Such a 
problem is beyond the scope of this study however. 

We took into account a nonzero electric field near the slab surface created by a double charge layer present at the 
surfaces of solids. It was found that, because of the presence of this electric field, the continuity equation for the 
axial current is fulfilled for the low-energy modes considered in this study without taking into account the linear in 
electric field anomalous term. In all fairness, the current understanding of the chiral anomaly in the problem of a 
finite-size Dirac semimetal remains incomplete. While it is understood that the chiral symmetry is explicitly broken 
by the surface effects, it is not completely clear whether the model implementation of the boundary layer used in the 
present study is sufficient to capture all relevant physics effects. 

One of the limitations of this study is the use of a special orientation of the magnetic field perpendicular to the slab 
surfaces. The main reason for this was that the analysis would become much more challenging in the case of a tilted 
magnetic field. From the physics viewpoint, the complications come from the need of squeezing the Landau orbits into 
a finite thickness of the slab. Nevertheless, here we could speculate that, in the case of a general orientation of the 
magnetic field, one may consider separately the perpendicular and parallel (with respect to the surface) components of 
the currents, while the results for the perpendicular components of the currents are expected to remain qualitatively 
the same as in the special case considered in this paper (except for the replacement B ^ B cosO, where 0 is an angle 
between the direction of the magnetic field and the normal to the slab surface). The parallel components of the 
currents are most likely to be the same as in an unbounded infinite system, but with B ^ B sin 9. While plausible, 
such a scenario clearly requires a further in-depth investigation, which is beyond the scope of this study. 
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Appendix A: Derivation of the Landau-level wave functions 

In this appendix we derive the Landau-level wave functions in the model described by the Hamiltonian in Eq. ©■ 
We look for the solutions of the Dirac equation, Htpir) = Etfj{r), in the form ip{r) = In this case, 

function satisfies the following equation; 

[{—vfPzJ^ — vpPy'l^ + ivpj^dx) 7 ° — vpeBxj^'y^ + 7717 °] cj){x) = Ecj){x). (Al) 

In the case of a constant magnetic field in the +z direction, assumed here, we choose the vector potential in the 
Landau gauge, A = (0, Bx, 0). (Clearly, all observables should be independent of a specific gauge choice.) Instead of 
the X coordinate, it is convenient to introduce a new dimensionless variable. 


(A2) 
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and rewrite Eq. dSB) as follows: 




- 7 = {PzY + {m/vp) - {E/vFh°) 

Vle^l 


HO = 0 , 


(A3) 


where sj_ = sign(eB). By making use of the same approach as in Ref. [h^ . let us introduce a set of linearly independent 
bispinors uf that satisfy the following relations: 


= Tut, 


[eH -VFP.H)ut = ssJvIpI - E^ u 

7 X _ ry2yp,p^ — TTl) 


± 

S 1 


■ = U' 


y/E'^ — Vppl — 

Then, by expressing the wave function as a linear combination of the two bispinors, 

M0 = Ki0ut + <^H0u;, 


(A4) 

(A5) 

(A 6 ) 


(A7) 


and using the relations in Eqs. (IA4I) through (IA 6 I) . we rewrite the equation for the wave function (IA3I) in the following 
form: 


T = 0 , [5| T + H] = 0 , (A 8 ) 

where k = y/E^ — VpPl — m? jf-h and = VF\/\e^\. Without loss of generality, let us choose s_l = sign(ei?) = +1. 
Then, the solutions of Eq. (IA 8 I) are given in terms of the parabolic cylinder functions [62| : 


$l(e) = D,y^{V20, niO = ^77,V2-i(v^e)- (A9) 

By requiring that the wave functions are finite at |^| —>■ 00 , one finds that k^/2 = n where n = 0,1,2... are 
nonnegative integers. In this special case, the parabolic cylinder functions Dn(p/^0 can be expressed in terms of the 
Hermite polynomials E[n{0'- En{V^0 = By making use of the definition of k, we also obtain the 

corresponding Landau-level energies: E^ = Vppl Tm^ T 2ne|. It should be noted that function is not finite at 

|(^| —>■ 00 when n = 0. However, by using the relation in Eq. (IA 6 I) . we can express through uj, i.e.. 


u 


+ 

S 


—my® — "fysy/vppl — E^ 


mTsy/El- VfpI ^_ 
V^n- UpPl - m? 


(AlO) 


Then, in the case of the LLL (n = 0), we see that the coefficient in the last equation vanishes if s = — 1. This means 
that, for the LLL, only one value of the spin, s = —1, is allowed. Eor the higher Landau levels, both spin projections, 
i.e., s = ±1, are allowed. The Landau-level wave functions are given by 


(l^siO = 


2 -r /2 


-HniO - n 


. m T Sy/m? T 2ne|^ e 




V 2 ^ 


Hn-liO 


where, according to Eqs. (IA4I) through (IA 6 I) . 



ut = 


mTsy/El -v^pI 


( -X4 


X2 

V 0 


u\pi - m2 

Einally, the explicit form of the Landau-level wave functions reads 


X4 = 


- ^pPi 

En - VFPz 


X2- 


^(r)„=o = 

= e^^^^+^yyy[Cx<j}n-s=+i{0 + C2H-,s=-im- 


(All) 


(A12) 


(A13) 

(AM) 
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where 


m 


n>0,s—-\-l 


= Cl 


( m+^m^+2nel^^ t^\ \ 

* En-vpp, 

y^{i) 


m+^m?+2n<il 


m 


n>0,s= —1 


= c. 


y/m^+2nel ^^ , , 

\ E^-vfP:. -'n'.s; 

/ „■ V™^+2’^'=i m-^m?+2nt\^ 

ynio 


(l){^)n=0,s=-l — Cq 


V 
/ 

V 


. m-^m^+2ne| 

\/m^+2n(.l ... 

En—VFVz 

0 \ 

0 

~ Eo-vpp^y^^^) y 


/ 


(A15) 


(A16) 


(A17) 


_^2 ! <2 

Here l^n(C) = f (^) are the harmonic oscillator wave functions. By making use of Eqs. (jA13l) through (IA17I) . 

Y 2'^n\yf^ 

it is straightforward to obtain the wave functions ([3]) and (|4]) in the main text. 


Appendix B: Matching solutions at the semimetal boundaries 

In this appendix, we derive the explicit expressions for the wave functions in the semimetal with a slab geometry 
by matching the general solutions, presented in Eqs. and ([n|) in the main text, with the corresponding vacuum 
solutions. This is done by making use of the Bogolyubov bag model, in which the wave functions outside the semimetal 
satisfy the same type Dirac equation, but with a large vacuum band gap M. 


1. Matching the wave functions of the n = 0 modes 


Let us first consider the n = 0 modes. In order to obtain normalizable solutions outside the semimetal, in the regions 
with z > a and z < —a, in Eq. ([3]) we will replace Pz with ip'^ and —respectively. With such a parametrization 
of the vacuum solutions, the energy of the n = 0 modes is given by A = ±-\/M^ — {vfp'zY'- In the limit Af —>■ oo, 
the corresponding solutions must have the same (finite) energy as the n = 0 modes in the bulk. This implies that 
vfp'z ~ M — (m^ + Vppi)/{2M). Therefore, to leading order in the inverse powers of M, the corresponding wave 
functions in the vacuum are given by 


V'(l*)n=0.2>a 


0,2<—a 




^-M(z-a)/vF 





M {z-\-a)/vp 



(Bi) 


(B2) 


By making use of Eqs. ®, (EH), and (IB2I) . we find that the matching conditions on the two sides of the semimetal 
slab, 


^'siab(a:, y,a)„=o = 
4'siab(a:, y,-a)„=o = 

lead to the following equation for pz'. 

vfPz cos {2apz) + 


C)n—0,z>a: 

(B3) 

'ip{x,y,-a)n=o,z<-a, 

(B4) 

msin {2apz) = 0, 

(B5) 
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where Pz 0. Note that = 0 is not allowed because it corresponds to a trivial solution for the wave function. 
Satisfying the boundary conditions also fixes the value of constant Cq in the superposition of two counterpropagating 
waves in the bulk solution (l5|); 


^ 2 ivra^ Eo{Eo+VfPz) 

Co =—e Co—z -^ 

m [m — ivfVz) 


(B6) 


By taking this into account, we derive the following final expression for the wave function inside the semimetal: 


where 


^siab(r)„=o = 


( 0 

2giap^ VFPi cos[p;^(z-a)]-{m+iEo) sin[pa(z-a)] 
im-\-VFPz —Eq 
0 

, VFPz cos\pz(z-a)]-(m-iEo) sin[p^(z-a)] 

\ im+VFp^-Ea / 


ICoP 


f-L [ry + vfPz{vfPz - -Eo)] 
Avf [2a(m^ + VpPz) + mvF] 


is obtained from the condition of the wave function normalization. 


(B7) 


(B8) 


2. Matching the wave functions of n > 0 modes 


Let us now consider the modes with n > 0 outside the slab. As in the case of the n = 0 mode, in order to obtain 
normalizable solutions in the z > a and z < —a regions, in Eq. o we replace Pz with ip'^ and —ip'z, respectively. 
The energies of such solutions are given by + 2ne^ — [vfp'zY- In the limit M —>• oo, these should 

coincide with the corresponding expressions for the Landau-level energies in the bulk. This is satisfied if we choose 
vfp'z ~ Af — (m^ + Vpp1)/{2M) in the vacuum solutions. By taking this into account, we derive the following vacuum 
wave function in the z > a region: 




I — p-p'z(^-°-)+VyV 

<n,z>a — ^ 


M {C'p + C'_) 


Yn-liO 

\/2^ 




VoJ 



C'+Y^iO 

0 

^C'_Y^{0 




+ 0 


(B9) 


where we kept terms up to subleading order in the inverse powers of M. From the normalization of the wave function 
(IB9I) . we find that 




2 M 



Therefore, the final expression for the wave function in the vacuum region z > a is given by 


V'(r)n,.>a = 

C[Y^-iiy 

o' 

1 


fn 









(BIO) 


(Bll) 


The wave function in the other vacuum region, z < —a, can be obtained in a similar way. The final result reads 


yvyz<-a = 


fol 


0 




\ 0 } 





(B12) 


The matching conditions at the boundary of the semimetal are similar to those in Eqs. (IB3I) and (IB4I) . i.e., 

4'siab(a;,2/, a)„>o = i’ix,y,a)n>o,z>a, 

^'siab(a;,2/,-a)„>o = i’ix,y,-a)n>o,z<-a- 


(B13) 

(B14) 




























16 


By substituting the wave functions from Eqs. (ID), (IBllI) . and (IB12I) into these matching conditions, we find that Pz 
should satisfy the same spectral equation (IB5I) as in the case of the n = 0 modes. By taking into account, however, 
that the rank of the system of Eqs. (IB13I1 and (IB14I1 is 2 units less than the dimension of the system, we obtain the 
following two linearly independent solutions for each higher Landau level (n > 0) inside the slab: 




^Sb(r)n>o = 


VFVz cos\pz(z+a)] + (m-iEn) sin[p^(z+a)] y (e\ 
m-\-i{vppz—En) 

.y/2n£|_6“*^*+£;„)] sin[p2(z+a)] y 

m^-\-ivppzm+2ne^ 

n„-iavz 'VFPz cos\pz{z+a)] + {rn+iEn) sinfp^(z+a)] 

\ m+i{vFP.-E„) / 

y^2rae|^y^m2+2ne|^ ” ^ 

r) e~’-^^‘^(vFPz+E„) sin[p^(z+a)] 

^ y/m^+2uel 

* -y/m^+ 2 ne|^[pa cosfpz {z-\-CL)]-\-{m—iEn) sin[pa (2+a)]] 


(B15) 


2i- 


V 


zj2ne\{vFPz -E-rP) 

-iav, (^>FP»+gn)sin[p;^(z+a)] y 
^Tn^+2nel 


-Yn-liO 


(BIG) 


where the normalization constants are given by 


|C+P 

|C_p 


clPz + nej + vfPz (vfPz - En)] 

Avf [2apzE^ + mvFPz — nt\ sin (dap^)] ’ 

_ ne\pz {vfPz - Enf _ 

+ 2ne\) [2apzE‘^ + itivfPz — ne\ sin (dap^)] 


(B17) 

(B18) 


3. Matching of the LLL wave functions in the case of finite vacuum gap 

Let us now consider the LLL in the case of large but finite vacuum gap M. In order to derive the tails of the LLL 
wave functions in the regions z > a and z < —a outside the slab, we use Eq. with pz replaced with arid —ip'z^ 
respectively, i.e., 


'4^{^'}n—0,z>a 


'?/^(r)n—0,2<—a 




M 


/ 




V 


V Eo—ivpp'z / 

0 

M 


Eo+ivpp'z 


(B19) 


(B20) 


Here p^ = Vf^V ~ • Enforcing the matching conditions at the boundaries of the slab, see Eqs. (|B3p and (IBdp . 

we find that the wave vector pz should satisfy the following spectral equation: 


-vfPz\/M'^ - El cos {2apz) + (Eg - mM) sin (2ap^) _ ^ 

M “ 


(B21) 
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[As expected, in the limit M —>■ oo, this equation reduces to Eq. (IB5I) .] The corresponding LLL wave functions inside 
the semimetal and outside the slab are given by 


/ 


^siab(r)„=o = AoYoiOe^P-y 


0 


^ip^{z-a) _ -ip^{z-a) 


-a) (-E'O+i'fPz) 


0 


m£'o(m—M) 


Y vfPz — Eq 


oW. 


_ ^-ivJz-a) {Eo-VFP.)\El-mM-ivpp,^M^-El 


\ 


/ 


, (B22) 


V'(r)„=o.z>a = 


VfPz 


mEQ{m—M) 

m{m - M) + {vfPz + Eq){vfPz - i\/M^ - E'^) 


mE^irn — M) 


; \ 

0 

M 


(B23) 


V'(r)«=o,.<-a = 


Egm 

VfPz cos {2apz) + \/— E^ sin {2apz) 


{Eg - VfPz) 

m{Eg - iy/M'^ - Eq) - M{vfPz + Eg) 


/ 




0 
1 
0 

\ EoFi^M'^-El ) 


(B24) 


where the overall constant Ag is obtained from the normalization condition. Its explicit expression is given by 

€Lim - M)Eo{Eo - VFPz) 


l^ol^ = 


Avf ( 2aEl[m - M) - mvF^/M^^g- -^= 


(B25) 
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